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We study the robustness, against the leakage of bosons, of wave functions of interacting many 
bosons confined in a finite box, by deriving and analyzing a general equation of motion for the 
reduced density operator. We identify a robust wave function that remains a pure state, whereas 
other wave functions, such as the Bogoliubov's ground state and the ground state with a fixed 
number of bosons, evolve into mixed states. Although these states all have the off-diagonal long- 
range order, and almost the same energy densities, we argue that only the robust state is realized 
as a macroscopic quantum state. 



When a quantum system is subject to perturbations from its environment, most wave functions decohere, and only 
an exceptional wave function(s) remains pure. For quantum systems with a single degree of freedom (/ = 1), this 
robust wave function is a coherent state 1^,^]. For example, when a coherent state \a) of single-mode photons passes 
through an absorptive medium, the final state is also a coherent state \oc'), which was attenuated (\a'\ < \a\) by the 
absorption Q|. It was also argued for a / = 1 system that coherent states produce the least entropy in the environment, 
thus being stable [Q. Since these conclusions are based on analyses on / — 1 systems, a natural question is: Are they 
applicable to macroscopic systems, i.e., to / 3> 1 interacting systems? Moreover, we must identify which coherent 
states are robust, because there are many choices of the coordinate (among many degrees of freedom) by which a 
coherent state is defined. Furthermore, for massive bosons the superselection rule (SSR) forbids superpositions of 
states with different numbers of bosons. Hence, we must clarify how coherent states can be compatible with the 
SSR. The purpose of this Letter is to answer these questions for condensates of interacting many bosons, which (or, 
equivalents of which) are observed in many physical systems such as liquid He ||] , quantum Hall systems Q , excitons 
||, and trapped atoms |6|. We also discuss the symmetry breaking in view of the robustness. 

We consider many bosons which interact with each other repulsively. We assume that the bosons are confined in a 
large, but finite box of volume V, which is placed in a huge room of volume V 3> V, which we call the environment. 
Suppose that the potential of the walls of the box is not high enough, so that the box and the environment exchange 
bosons via tunneling processes at a small rate (flux) J. Let t eq denote the time scale after which the total system, 
the box plus the environment, reaches the equilibrium state. We are not interested in the time region t > t cq because 
the equilibrium state is just a uniform state that is determined solely by the initial state of the environment (because 
V V). We therefore examine the transient region for which t <C t eq , in order to discuss the robustness of an initial 
state |0(O)), which is prepared at t = 0, of the box. Depending on the choice of the initial state |</>e(0))e of the 
environment, the box state may be affected either drastically or moderately. For example, a moderate situation is 
that |0e(O))e has the same density n of bosons as |</>(0)). In such a case, n of the box will be kept constant for all 
t. To discuss the robustness, however, we consider the severest situation where the environment is initially in the 
vacuum state |0)e of bosons, so that bosons escape from the box continuously. If a box state is robust in this severest 
case, it would also be robust in other cases. Hence, the total wave function at t = is | < &(0)) to ta.i = 1 0(0)) ® |0)e- We 
decompose (the r dependence of) the boson field that is defined on V + V as V>totai(r) = V>(r) + V'E(r). Here, ^(r) 
localizes in the box, whereas the low-energy component of V>e(f) localizes in the environment 0. Accordingly, the 
Hamiltonian of the total system is decomposed as H tota \ = H + He + i?sE- Here, H (He) is a function of y> ($e) only, 
describing interacting bosons in the box (environment). On the other hand, Hse includes both ip and ipE, describing 
the ip-'tpE interaction. If the leakage flux J is small, the probability of finding two or more bosons simultaneously in 
a wall of the box is negligible, and thus the dominant term of Hse takes the following form: 



Here, w(r) represents the shapes of the walls (w ~ 1 in the walls, w = in other regions), whose potential height is 
characterized by a parameter A. Details of w are irrelevant because they are all absorbed in the value of j, Eq. (||). 
In the time region of interest (t -C t eq ), the ipE-tpE interaction should be unimportant because n of the environment 
remains zero. On the other hand, we must treat the tp-ip interaction appropriately. For this purpose, we use the 
decomposition formula for ip &[Hj] : 
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$ = S + ^', (2) 

where S is an operator satisfying 

&\N,G)=VN£\N-1,G), (3) 



where \N, G) denotes the ground state that has exactly N bosons which we call the number state of interacting 

bosons (NSIB), and 



t = (N-l,G\i/>\N,G)/VN (4) 



is a hallmark of the condensation: \A/V£ = 0(1) for condensed states, whereas y/N£ = 0(1/ y/V) for normal 
states ||,|9|,|l2|]. We here consider the condensed states. Since ip alters A" exactly by 1, Eq. ([!]) means that 
(N - AN, G\i>\N, G) = VNZ 5 A n,i for all AN such that \AN\ < N. It then follows from Eqs. (§) and (§) that 

(N — AN, G\tp'\N, G) = (for \AN\ < N). (5) 

Namely, ip' transforms \N, G) into excited states. 

For weakly-interacting bosons, the explicit forms of the NSIB were given in Refs. §,0, and that of a was 
given in ||. Because of the boson-boson interaction, they are rather complicated functions of bare operators 
a k : |JV,G) = (l/y/Nl.y d (a f o ) N \0) and 5 = n Q / nV e id a e~ ld . Here, n Q = (N - J d 3 r^<tp')/V, G = 

(—i/2nV)a\a\ ) ^ ^ n y g a q a_ q + H.c, ip is an arbitrary phase, and y q is given in Ref. Using these expressions, we 
can show that HQ] 

[S,St], [S, ip'], [S,^t] =0(1/V). (6) 

Lifshitz and Pitaevskii (LP) |8| claimed that Eq. (^) is applicable even when the interaction is stronger. Their 
discussion is somewhat controversial because LP started from, instead of Eq. (||), the assumption that S could be 
defined by 3|JV, v) = "E^N — l,v), where \N , v) denotes any eigenstate that has exactly N bosons. However, we note 
that for weakly- interacting bosons we have not used this assumption in the derivation of Eq. (||) . We thus expect that 
Eq. (^|) also holds for bosons with stronger interaction, even if LP's assumption was too strong. If this is the case, the 
following results are applicable not only to weakly interacting bosons but also to bosons with stronger interaction, 
because the results will be derived only from Eqs. (0)-(||). 

Since we are studying the robustness against weak perturbations, we assume that A is small, so that J is very small. 
In this case, we have to consider transitions only among \N, G)'s with different A^'s (i.e., we can neglect transitions 
to excited states). Hence, the reduced density operator p can be generally written as 

P(t) = Y J PNM(t)\N,G)(M,G\. (7) 

N,M 

It seems almost obvious that quantum coherence between \N, G) and \M, G) with large \N — M\ would be destroyed 
by the interaction with the environment. We therefore study the most interesting case where pnm is localized in the 
N-M plane in such a way that yj (5N 2 ) <C (N). If this relation is satisfied at t = 0, it is also satisfied for all t -c t cq . 
We also assume that V is large enough, so that the boson density in the environment (= [(A"(0)) — (N(t))]/V) is 
negligibly small for all t <§; f eq - Under these conditions, we can calculate the time evolution of pNu(t), using Eqs. 
(D-(|, as @U 

p NM (t + At) = e -i(N-MMn{t))At/h 

x [ P NM(t)(l — (N + M)j(n(t))At/2) (8) 



+ p N +i,M+i(tW(N + 1)(M + l)j(n(t))At\ + 0(\ A ), 

for a finite time interval At that satisfies Ti/E c < At < 1/ (N)j(n), where E c is the energy scale over which the matrix 
elements of Hse are non-negligible. Here, n = (N)/V, and p(n) (> for a condensate of interacting bosons |||],D) 
denotes the chemical potential of bosons in the box. Furthermore, 

j(n) = K 2 ^\\\^D(»(n)), (9) 
n n V 
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where D(p) is the density of states per unit volume of the environment at energy p, v is the total volume of the walls 
of the box, and K is a constant of order unity. Both p and j depend on (N) through n, but this dependence is very 
weak because a change of (N) by 1 only causes the change of n by l/V. Note that our basic equation (||) has only 
two parameters, p and j. Namely, all model-dependent parameters (details of H , and H^e) are absorbed in these 
two parameters. Therefore, the following results are general and model independent. 

Using Eq. (^), we first calculate the time evolutions of the expectation value (N) and the fluctuation (SN 2 ) of the 
number N of bosons in the box. We find 

±{N) = -j(n)(N). (10) 

Hence, (N) decreases gradually because of the leakage flux J — j(n)(N). For (8N 2 ), on the other hand, we find 

±F = j(n)[l-F], (11) 

where F = (6N 2 )/(N) is the "Fano factor" §. It is seen that a robust state must have F = 1, whereas any states 
with F ^ 1 are fragile in the sense that their F evolves with time, approaching unity. For example, the ground-state 
wave function in the Bogoliubov approximation, |Bog, G), has F > 1 ||. Hence, it is fragile. The ground-state wave 
function with a fixed number of bosons, \N, G), is also fragile because F = 0. 

Since the evaluation of F is easy, F is a convenient tool for the investigation of the robustness. However, since 
F is only related to the diagonal elements of p, it does not distinguish between pure and mixed states. Therefore, 
we now solve the basic equation (^) for various initial states to investigate the robustness of the wave functions in 
more detail. When the initial state is a pure state of the NSIB, i.e., p(0) = \N, G)(N, G|, then pnm after a short 
interval At is evaluated as pNN{At) = (1 — Nj(n(t))At), pN-i,N-i{At) = Nj(n(t))At, and other elements are zero. 
Therefore, p becomes a classical mixture of \N, G) and \N — 1, G) at t = At, in consistency with the above result that 
states with F = are fragile. By evaluating the evolution at later times, we find that p evolves toward a Poissonian 
mixture of \N, G)'s |0|, consistent with F — ► 1. In a similar manner, we can show that the pure state of Bogoliubov's 
ground state p(0) = |Bog, G)(Bog, G|, which has F > 1, also evolves into a mixed state. We can also show that the 
number-phase squeezed state of interacting bosons (NPIB), which was found in Ref. Q as a number-phase minimum 
uncertainty state with < F < 1, also evolves into a mixed state. These examples show that F is indeed a simple 
measure of the robustness: A pure state with F ^= 1 is unlikely to remain pure. Note, however, that a pure state with 
F = 1 is not necessarily robust. For example, we can show that the coherent state of free bosons (CSFB) evolves into 
a mixed state, although it has F = 1. Hence, F = 1 is only a necessary condition for the robustness. 

Among many states with F = 1, we have successfully found a very special state that is robust in the sense that it 
remains pure when it is weakly perturbed by the environment. The state is given by 

p(t) = \a(t),G){a(t),G\. (12) 

Here, a(t) is a time-dependent complex number given by 



a(t) = e^./(NW), (13) 

where (N(t)) is the solution of Eq. (|io|), and 

<p(t) = <p(0) - l - \ p(n(t))dt. (14) 



h Jo 

Here, the initial phase <^(0) is arbitrary, and n[t) = (N(t))/V. Furthermore, 

00 M 

|a,G)^e-H 2 / 2 ^ ^=|Af,G), (15) 

Af=0 V M - 

which we call the coherent state of interacting bosons (CSIB). It has the same form as the CSFB except that 
\M, G) is the NSIB. Because of this difference, simple relations for the CSFB do not hold for the CSIB. For ex- 
ample, (a,G\^}\a,G) =/= a/^/V, and, moreover, \a, G) is not an eigenstate of Nevertheless, {a, G\N\a, G) = 
(a,G\SN 2 \a,G) = \a\ 2 , hence F = 1 exactly, as in the case of CSFB. Since the NSIB has a complicated wave func- 
tion, so does the CSIB. [For weakly-interacting bosons, its explicit form was given in Ref. 0.] Although complicated, 
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the wave function of the CSIB is robust against weak perturbations from the environment: It keeps the same form, 
whose parameter a(t) evolves slowly (except for the phase rotation), and remains a pure state, in contrast to other 
wave functions which soon evolve into mixed states. In fact, Eqs. (|12|)-([15|) yield 



e i(N-M){ip—nAt/h) 



X 



VMW. 

, . , N + M , N + AI 

(JV(t))->-(l-iAt)-=-, (17) 



which indeed satisfy Eq. (g). 

We now discuss the compatibility with the SSR, which might raise the objection that the CSIB would not be 
realized because superpositions between states with different values of N are forbidden for massive bosons. To show 
that this intuitive objection is wrong, it is sufficient to give one counterexample. Suppose that there is another box, 
which also contains condensed bosons, in the same room. The total system consists of two boxes and the environment. 
According to the SSR, the wave function of the total system |$) total should be a superposition of states that have the 
same number of bosons, iVtotal = N + N' + N-g = fixed, where N' denotes the number of bosons in the second box. 
Consider the following state, which satisfies this constriction; 

N,N',e 

x \N, G) ® \N', G)' <8 l^totai -N-N', £) E - (18) 

Here, a = |a|e lv , a' = |a'|e iV , and Cg is a complex number, where I is a quantum number labeling states of the 
environment \M,£)e which has M bosons. Regarding the phases cp and <p', only the relative value ip — ip' = 9 has a 
physical meaning. We thus take ip' = henceforth. Equation (|l8|) yields the reduced density operator of the first box 
as P = J2n e ~ |a|3 (\ a \ 2N /iV!)|iV, G><iV, Gj. It is easy to show that this is identical to 

— ||a|e",G><|a|e* G|. (19) 

Although this p represents a mixed state of CSIB's, we note that it does not contain the maximum information on the 
state in the box, whereas the best density operator should have the maximum information. The lacking information 
is that the phase relative to the condensate in the second box is ip. Hence, the maximum information is Eq. jl9| ) 
with the restriction 9 = ip. This combined information is concisely expressed as p = ||a|e lv , G)(|a|e lv , G|, which 
agrees with Eq. (|lj). Namely, Eq. ( |l2] ) is better than Eq. (19) because it contains more information. This example 



demonstrates that Eq. (|12j) can be compatible with the SSR in realistic cases where the box exchanges bosons with 
the environment. Only in the limiting case where the box is completely closed, should the SSR be crucial, and the 
NSIB would be realized if the temperature T — > H] . 

We have established that the CSIB is a robust pure state of interacting many bosons. We finally discuss its 
implications. The robustness of the present work should not be confused with the "stiffness of macroscopic wave 
functions" (^,^1, which only referes to the stability of an order parameter in a mean field approximation. For example, 
| Bog, G) has the stiffness HQ , whereas it is fragile as we have shown. The robustness is generalization of the robustness 
of coherent states of / = 1 systems jlj|2|. It is thus natural to expect for / ^> 1 systems that some coherent state 
would be robust. However, it was not known which coherent state is robust: there are many choices of the coordinate 
by which a coherent state is defined. Since Eq. (|^) yields 

(E/0\a,G)=a\a,G), (20) 

the present work has revealed that the robust coherent state is the one defined by S/£. In this sense, S + 5^ is 
the "natural coordinate" of interacting many bosons. The condensation of bosons are often characterized by the 
off-diagonal long-range order (ODLRO) that is defined by (^t(r)^(r')) = finite for |r- r'| - V 1 / 3 @. Using Eqs. 
(|)-(|) we can show that the CSIB, NSIB, NPIB, and the Bogoliubov's ground state all have the ODLRO. Hence, 
the present work has revealed that the ODLRO does not necessarily imply the robustness. Furthermore, all of these 
states have almost the same energy densities, i.e., the differences of (H)/V are only 0(1/V) for the same value of 
(N) Q. For example, if we let En,g be the eigenenergy of the NSIB, H\N,G) = En,g\N,G), we can then easily 
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show from Eq. tfljfy, neglecting terms of 0(l/V), that (a,G\H\a,G)/V = E H 2^ G /V = E {N):G /V. Therefore, the 
robustness of the CSIB is not due to an energy difference, but to natures of wave functions. Since interactions with 
the environment are finite in most physical systems, we argue that only the robust state, CSIB, should be realized 
as a macroscopic pure state. Since the (relative) phase of the CSIB is almost definite ||, the global gauge symmetry 
is then broken. Although V is finite, we are thus led to the symmetry breaking by considering the robustness. This 
suggests that quantum phase transitions may have more profound origins than singularities that are developed as 
V — ► oo. A conventional trick to get symmetry breaking states for boson condensates is to introduce a symmetry 
breaking field r\, which couples to ip as H v — J d 3 r(ri*tf; + r]tf;' i: ). However, 77 is usually considered as an unphysical field 
P, [l6|]I^] , and it was sometimes argued that symmetry breaking states were meaningless because they look against the 
SSR Jl7[ ] . In contrast, the present work gives a physical reasoning for the symmetry breaking, assuming only physical 
interactions, and shows the compatibility with the SSR. 
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